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O ■ ABSTRACT 

o 
o 

I We investigate the effect of a hybrid electron population, consisting of both thermal 

^Lf, and non-thermal particles, on the synchrotron spectrum, image size, and image shape 

I of a hot accretion flow onto a supermassive black hole. We find two universal features 

in the emitted synchrotron spectrum: (i) a prominent shoulder at low 10^^ Hz) 
frequencies that is weakly dependent on the shape of the electron energy distribution, 
and (ii) an extended tail of emission at high 10^^ Hz) frequencies whose spectral 
slope depends on the slope of the power-law energy distribution of the electrons. In 



> 

Q> , the low-frequency shoulder, the luminosity can be up to two orders of magnitude 



o 
o 



greater than with a purely thermal plasma even if only a small fraction (< 1%) of 



I the steady-state electron energy is in the non-thermal electrons. We apply the hybrid 



model to the Galactic center source, Sgr A*. The observed radio and IR spectra imply 
that at most 1% of the steady-state electron energy is present in a power-law tail in 





that such a hybrid distribution can be sustained in the flow because thermalization via 



Q^' this source. This corresponds to no more than 10% of the electron energy injected 

Q . into the non-thermal electrons and hence 90% into the thermal electrons. We show 

I Coulomb collisions and synchrotron self-absorption are both inefficient. 

The presence of non-thermal electrons enlarges the size of the radio image at low 
^ frequencies and alters the frequency dependence of the brightness temperature. A 

^ • purely thermal electron distributions produces a sharp-edged image while a hybrid 



distribution causes strong limb brightening. These effects can be seen up to frequencies 
~ 10^^ Hz and are accessible to radio interferometers. 
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1. Introduction 

The mechanisms of particle heating and acceleration, and the emission spectra from the 
resulting particle energy distributions, are of great importance in the theory of collisionless hot 
accretion flows onto compact objects. Discussions in the literature have focused on the physics 
of electron heating and acceleration (Begelman & Chiueh 1988; Bisncmityi-Kogan & Lovelace 
1997; Quatacrt &: Gruzinov 1999; Gruzinov & Quataert 1999, Mcdvcdev 2000), the efficiency of 
particle thcrmalization (Ghisellini, Guilbert, & Svensson 1988; Ghisellini, Haardt, & Fabian 1993; 
Mahadevan & Quataert 1997; Ghisellini, Haardt, &; Svensson 1998; Nayakshin & Melia 1998), 
and the generation of hybrid thermal-nonthermal electron energy distributions in these plasmas 
(see, e.g., Coppi 1999 and references therein). Despite a substantial amount of work, many issues 
remain unresolved, primarily because of our incomplete understanding of physical processes such 
as magnetic reconnection, MHD turbulence, and collective plasma modes. 

These questions are especially relevant for an optically-thin advection-dominated accretion 
flow (ADAF). An ADAF is an example of a hot, rarefied, magnetic plasmas with low radiative 
efficiency (Ichimaru 1977; Narayan &: Yi 1994, 1995b; Abramowicz et al. 1995; see Narayan, 
Mahadevan, &; Quataert 1998b and Kato, Fukue, &; Mineshige 1998 for reviews). A basic property 
of the nearly collisionless plasma in an ADAF is that the Coulomb coupling between the electrons 
and ions is weak so that energy transfer from the ions to the electrons is inefficient. In addition, it 
is commonly assumed that the viscously generated energy primarily heats the heavier species, the 
ions, and that the electrons retain a thermal distribution throughout the flow (e.g., Narayan &; Yi 
1995b; Mahadevan 1997). 

There are, however, processes such as MHD turbulence, pair production (e.g., through pion 
decay), and electron-proton coupling which can both heat the electrons and generate non-thermal 
distributions. Quataert and Gruzinov (1999; see also Gruzinov k. Quataert 1999) considered 
two processes specific to MHD turbulence that accelerate particles in magnetic collisionless 
plasmas: Landau damping by electric fields parallel to the local magnetic field and transit-time 
damping by time-varying magnetic fields. They found that the assumption of negligible electron 
heating/acceleration is valid only for weak magnetic fields, i.e., when the ratio of the gas pressure 
to total pressure /3adaf is larger than a critical value /3crit- The value of /3crit is very uncertain and 
is around 0.9. For /?adaf ^ 0.9, turbulence primarily accelerates the protons, while for stronger 
magnetic fields, the results are inconclusive. Shocks and pion decay can also lead to non-thermal 
electrons in the accretion flow. 

Several processes, such as Coulomb collisions and synchrotron self-absorption, can potentially 
lead to thermalization of particles in accretion flows (Mahadevan k, Quataert 1997; Ghisellini 
et al. 1998). Mahadevan & Quatacrt (1997) showed that Coulomb collisions are ineffective in 

thermalizing the electrons in an ADAF. However, they argued that, for sufficiently high mass- 
accretion rates, the electrons in the plasma can be thermalized by synchrotron self- absorption. 
Nayakshin and Melia (1998) showed that considerable deviations from a Maxwell-Boltzmann 
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distribution can be sustained in a plasma with low source compactness when Coulomb collisions, 
Comptonization, and pair processes are taken into account. 

In view of the difficulty of calculating the heating, cooling and thermalization of particles from 
first principles, many authors have restricted their models to either purely Maxwellian or purely 
non-thermal (extended power laws or monoenergetic) electron distributions. Only recently have 
there been attempts towards explaining spectra of accreting black holes with models including 
hybrid thermal/non-thermal distributions of electrons. Some models have physically motivated 
distributions, such as non-thermal electrons produced by decaying pions (Mahadevan 1998), 
whereas others invoke more ad hoc distributions to fit the data (e.g., Beckert & Duschl 1997; 
Falcke k. Biermann 1999). The reverse process of trying to constrain the energy distributions of 
particles in accreting plasmas by comparing models to data, however, faces issues of uniqueness 
which can be addressed only by a more comprehensive study of hybrid models. 

In this paper, we consider generalized electron distributions consisting of a dominant 
Maxwellian plus a small non-thermal power-law component of varying slope and energy content, 
and study the synchrotron emission from the resulting hybrid plasmas. We identify the 
characteristic signatures of the non-thermal electrons on the emitted radio synchrotron spectrum 
of an accretion flow and on its image as observed with a radio telescope. We also discuss to what 
extent the observed effects could be used to determine the details of the underlying non-thermal 
electron distribution. This work is relevant for interpreting observations of low luminosity AGNs 
such as the Galactic Center source, Sagittarius A* (Sgr A*). 

In §2 we review the basic properties of ADAFs, hybrid plasmas, synchrotron radiation, and 
radiative transfer. In §3 we present a series of models of supermassive black holes with low 
accretion rates. We apply these models to Sgr A* and derive constraints on the fraction of the 
electron energy that can be present in a non-thermal form. In §4, we study the correspondence 
between the energy distribution of the electrons and the resulting spectrum. In §5 we study 
the energetics of a hybrid flow and calculate the heating, cooling and thermalization rates of 
the non-thermal electrons. We summarize our conclusions in §6. We present in an Appendix 
approximate analytic expressions for the contribution of a non-thermal particle population to the 
synchrotron spectrum of an accretion flow. 



2. Formalism 

2.1. Advection-Dominated Flows 

We begin by reviewing some of the basic properties of the optically-thin branch of ADAFs. 
ADAFs are quasi-spherical, hot, magnetic accretion flows in which the accreting plasma is too 
rarefied to cool efficiently by radiative processes. The viscously dissipated energy is therefore 
advected into the black hole or other compact object at the center (see Narayan et al. 1998b 
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and Kato et al. 1998 for reviews). In the limit where the fraction of the viscous energy advected 
inward is independent of radius, a self-similar analytic solution for the thermodynamic quantities 
of the accreting gas can be obtained (Narayan & Yi 1994, 1995b). We make use of this solution 
in Appendix A. For the numerical calculations presented in the rest of the paper, we use more 
accurate global solutions to obtain the run of electron temperature and density with radius. These 
solutions are calculated by the methods described in Narayan, Kato, k, Honma (1997b), Chen, 
Abraniowicz, k. Lasota (1997), and Popham & Gammie (1998). The magnetic field strengths are 
obtained by assuming that the ratio of gas pressure to total pressure (sum of gas and magnetic 
pressure) is equal to a specified value /3adaf- 

Since the focus of this paper is on massive black holes in galactic nuclei with low mass-accretion 
rates, with specific applications to the black hole in our own Galactic nucleus, Sgr A*, we scale 
masses in units of 10® solar mass, i.e., M = melO^M©, and radii in units of the Schwarzschild 
radius, i.e., R = rRsch-, where 

2GM 

Rsch = — ^ = 2.95 X lO^^me cm. (1) 

We scale the mass accretion rate in units of 10~^-MEdd) i-e., M = ?fi_3lO~^MEdd! where the 
Eddington mass accretion rate is 

MEdd = ^ = 1.39 X 10^^ (^y'me g s-\ (2) 

In defining the Eddington rate, we assume a standard radiative efficiency of r/cff = 0.1. (This is 
purely for the purposes of the definition; the actual radiative efficiency can be very different from 
0.1). 

For the calculations presented here we use ADAF models in which the viscosity parameter is 
set to a = 0.1, the equipartition parameter to Padaf = 0.968, and the ratio of viscous electron 
heating to proton heating to 6 = 10"^. Note that /3adaf differs from the usual plasma parameter 
/?, which is the ratio of the gas pressure to the magnetic pressure; the value of /3adaf = 0.968 
assumed in our ADAF models corresponds to a plasma /? = 10 (Quataert & Narayan 1999). 



2.2. Hybrid Populations 

We assume that a large fraction of the electrons in the plasma are in a thermal distribution 
with temperature T and that the rest of the electrons are in a non-thermal distribution, usually 
with a power-law form. We denote the number density of electrons in the thermal population 
by iVth and in the non-thermal population by A^pi. We denote the emissivities of the thermal 
and power-law electron populations by jth and j^i respectively, and the corresponding absorption 
coefficients by ath a^nd api. 
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For the thermal electron population, we use the relativistic Maxwell-Boltzmann distribution 
given by 

nth(7) = iVth7'/3exp(-7/0e)/[^eK2(l/^e)], (3) 

where 7 is the electron Lorentz factor, P is the relativistic electron velocity, and 9e = kT/mef? is 
the dimensionless electron temperature. The modified Bessel function of second order K2{'i/0e) 
arises from the normalization of the Maxwellian. Similarly, for the non-thermal electron population 
we use a power-law distribution extending from 7 = 1 to infinity, 

ripi(7) = iVpi(p-l)7-^ (4) 



The number density of thermal electrons Nth is a function of the flow radius and is determined 
by the global ADAF solutions. We determine A'pi at each radius by assuming that the steady-state 
energy in the non-thermal distribution is equal to a fraction 77 of the energy in the thermal 
distribution, with 77 constant in radius. Although the calculations presented in §3 are all carried 
out with this assumption, generalizations to radially-dependent r/(r) as well as a discussion of 
the energetics of such a flow will be presented in §5. Note that we implicitly assume that the 
non-thermal electron population does not affect the dynamics or the thermal properties of the 
flow; this will be justifled in §5. 

The energy density of a MaxwcU-Boltzmann distribution of electrons at temperature 6^ was 
derived by Chandrasekhar (1939, eq. [236]) to be 

U = a (6'e) NthrUeC^Oe, (5) 

where 



(6) 



4i^2(l/^e) 

is a coefficient that varies from 3/2 for a non-relativistic electron gas to 3 for fully relativistic 
electrons, and are the modified Bessel functions of the nth order. For the present purposes, 
we use a simplified expression for a{9e) which has an error of less than 2% at all temperatures 
(Gammie Sz Popham 1998): 

6 + IbOe 

The number density of the non-thermal distribution is then 

Npi = ria{ee)ee{p-2)Nth. (8) 
This normalization of the power law population typically corresponds to 

^ ^ (0.1 - 10)7/, (9) 
depending on the electron temperature and the power law index. 
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The distributions considered above correspond to the steady state that results from the 
competition between heating/acceleration and cooling by radiation. We discuss in some detail in 
§5 the energy equations for the thermal and non-thermal electron populations. Here we simply 
note that the synchrotron cooling timescale of electrons moving with a Lorentz factor 7 scales as 
7~^, and hence electrons in the high energy tail of a power law distribution cool most rapidly. 
As a result, if electrons are injected with a power-law distribution with index s [n{'y) oc 7"''] and 
cool only by synchrotron emission, the synchrotron cooling causes the power law index of the 
steady state distribution to be p = s + 1 above a certain 75, called the cooling break, thus causing 
the steady state distribution to fall off more steeply at higher electron energies. Mahadevan & 
Quataert (1997) calculated 7^ in an ADAF by comparing the cooling timescale to the inflow 
timescale and found that at sufficiently high accretion rates, the break occurs at a very low Lorentz 
factor, 7(, ~ 1.5. The Lorentz factors of interest to us are invariably larger than 7;,. Therefore, the 
values of p we consider below in §3 and §4 arc always equal to s + 1, so that the injected energy 
distribution 7"'' is harder by one power of 7 than the steady state energy distribution, 7"*^. 

Corresponding to 77, we can define another quantity r/inj that measures the fraction of electron 
energy injected into a power law distribution. If s < 2, then rji^j can be significantly greater than rj. 
If wc assume a distribution from 7min = 1 to some 7max5 ^inj is greater than r/ by a factor 7inax- 
However, it is thought that the acceleration mechanisms typically encountered in astrophysics, 
such as shock acceleration or acceleration via MHD turbulence, inject energy into particles with 
s > 2 such that the steady state distribution has p > 3. In this case there is little dependence on 
7max, and 

5^ = (10) 

which is not very different from unity. Although we expect s > 2 and p > S for most systems, for 
completeness we consider models in the range 2 < p < 4. 



2.3. Synchrotron Emissivity 

The synchrotron emissivity of a relativistic electron moving with a Lorentz factor 7 in a 
magnetic field of strength B is given by (Rybicki & Lightman 1979) 

juh,0) = ^'ybsmeFix). (11) 

Here, Uf, = eB /2TrmeC is the non-rclativistic cyclotron frequency, 6 is the angle between the 
direction of the magnetic field and the velocity of the electron, and 

00 

F{x) = X J K,/2{t)dt, (12) 

X 

with the modified Bessel function of order 5/3, x = v/vc, and Uc = ^j^UbSmO. For a thermal 
distribution of electrons, the total emissivity for a given angle 6 is obtained by integrating equation 



-7- 



3vM0) = r- ,^ -— , (13) 



(11) over the Maxwellian distribution (Pacholczyk 1970) 
where 

" 3^ ^^'^ 

and 

1 °° 

/(xm) = / z"^ e-x.p{-z)F{xM/z'^)dz. (15) 



The hmiting behaviour of I{xm) for small and large xm was derived by Pacholczyk (1970) and 
Petrosian (1981), respectively. Mahadevan, Narayan, &: Yi (1996, hereafter MNY96) integrated 
equation ([l3|) over the angle 9 for an isotropic distribution of electron velocities and provided a 
fitting function for the emissivity in the ultrarelativistic to mildly relativistic regimes: 



with M{xm) given by 



, , 4.0505 a / 0.40 h 0.5316 c \ , 1/3, , , 

M{xm) = Tu- 1 + ^TT + -^JT- exp(-1.8896xi/'). (17) 



1/6 I ' 1/4 ' 1/2 

The best fit values of the coefficients a, 6, and c for different temperatures are given in MNY96. 
The coefficients tend to unity in the ultrarelativistic limit. Finally, the synchrotron absorption 
coefficient ath is related to the emissivity via Kirchoff 's law, 

au,t\i= iu,t\i/By{T), (18) 

where Bjj{T) is the black body source function. 

For the total emissivity of electrons in a power law distribution, we use the expression given 
in Rybicki & Lightman (1979) and average over angles, 

i-,pi = Cl,r^-^a{ee)eeiy, iy-J , (19) 



where A^pi is defined in terms of A'th as above and 



, _ (P-i)(P-2) r(| + l|)r(|-3^)r(| + |) 

p'" 4 ip + i) ni + i) ■ ^'"^ 

The corresponding absorption coefficient is 

au,p\ = C^in ^ a{ee)Be\ — \ V (21) 

with 

/3^3P/2 (p - i){p - 2) r(^)r( ^)r(g±g) 

me r 



4 
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2.4. Radiative Transfer and Numerical Methods 

The equation of radiative transfer for a time- independent, sphericahy symmetric flow is (e.g., 
Mihalas 1978) 

fi{d/dr) + r"^(l - fi^){d/dfi)]I{r, fx, v) = j{r, v) - a{r, v)I{r, /i, v), (23) 

where = cos 6 = (dz/ds) is the cosine of the angle between the ray and the radial direction, r is 
the radial coordinate, is the frequency, and j and a are the emission and absorption coefficients 
defined above. One can simplify this equation by taking plane parallel rays of varying impact 
parameters (perpendicular distances of rays to the central line of sight) through the flow and 
solving the equation along these rays (Mihalas 1978). The equation then becomes 

f)T± 

±^=ju-ajt, (24) 

where now s is the line element along the ray and the coefficients +1 and —1 correspond to 
radiation coming towards and going away from an external observer, respectively. In our problem, 
3u = jth + jpi and = ath + Opi- Rewriting the equation in terms of the source function 
Sy = ju/oiy and optical depth r(s) = — ads' , where Sout is the point of intersection of the ray 



with the outer boundary of the flow, equation (24) becomes 



±^ = Iu-S,, (25) 

OT 



where the combined source function is 



r, 3v Jth + ipl 'S'th Sp\ , , 

Su = — = ■ = — -, + — — —, • (26) 

Oiu ath + ttpi 1 + Opi/ath 1 + ath/api 



In the numerical calculations reported below, we integrate equation ( p5| ) using the formal 
solution and the appropriate boundary conditions for a non-illuminated atmosphere (see Mihalas 
1978). We carry out the integral along rays with impact parameters up to ~ 2000 Schwarzschild 
radii, beyond which the temperature of the electrons becomes too low for signiflcant synchrotron 
emission. Because of the very steep dependence of the synchrotron emissivity on photon frequency, 
magnetic field strength, and electron temperature and density, we use an adaptive step size for the 
radiative transfer integral. The total flux is obtained by integrating over all impact parameters. 
We validated the implementation of our numerical algorithm by comparing its output to analytic 
solutions of the radiative transfer equation in uniform media. The discrepancy between the 
numerical and analytic solutions was <^ 0.5% for the cases considered. 

In all model ADAF spectra published so far, the transport of radiation has been calculated 
using an approximate method based on concentric shells (Narayan, Barret, & McClintock 1997a). 
In Figure 1, we compare, for a typical flow, the exact spectrum obtained by formally solving the 
above transfer equation to that obtained with the previous approximate method. We see that 
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there is a fairly good agreement between the two methods, but with some differences. The most 
prominent difference is a shifting of the peak toward higher frequencies in the exact calculation, 
as well as some broadening. This is probably due to the poorer resolution of the concentric shell 
method which becomes a limiting factor close to the black hole. There is also a slight offset at 
lower frequencies, probably again due to poor resolution. Note that there is a second peak at 
high frequencies 10^^ Hz) in the spectrum calculated with the approximate method. This 
is caused by the inverse Compton scattering of soft photons, a process which is not included 
in the radiative transfer code described here. The remaining features of the two spectra are 
quantitatively consistent. 




log (Frequency/Hz) 

Fig. 1. — The solid line shows a typical spectrum calculated with the radiative transfer method employed in this 

paper. The model parameters are those used for Sgr A*: mg — 2.5, m_3 = 0.1, and other ADAF parameters as 
specified in the text. The dashed line is an approximate spectrum of the same model, calculated by a simplified 
method described in Narayan et al. (1997a). The secondary rise at u > 10^^ Hz is due to Compton scattering, which 
is not included in the exact calculation. 



3. Numerical Results 
3.1. Parameter Study 

In this section, we study the effects of an extended power-law electron distribution on the 

synchrotron emission spectra of ADAFs. Figure 2 shows the various components of the spectrum 
of a typical hybrid model, with p = 3.5, i] = 1%, me = 2.5, and m_3 = 1. Compared to the 
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spectrum of a pure thermal model (dashed line), we see two primary effects due to the power-law 

electrons: (i) there is a prominent shoulder of optically thick emission at low frequencies and (ii) 
there is an extended power-law tail of optically thin emission at high frequencies. In between these 
two features there is a region of the spectrum where the thermal peak dominates. These features 
were first identified by Mahadevan (1998) for a specific model. We find that they are universal for 
hybrid models. 
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Fig. 2. — Regions of a hybrid synchrotron spectrum labeled according to the dominant emitting and absorbing 
electron population and the resulting spectral shape. The different segments and the transitions are present in all 
hybrid synchrotron spectra. 




To understand the results, we note that the spectra of hybrid populations are determined 
by a competition between jpi and jth and between api and a^^, each having different radial and 
frequency dependences. Since the emission at different frequencies arises from different radii in the 
flow, the relative importance of api to ath and jpi to jth at that particular frequency and radius 
determines the local behaviour of the spectrum. In the less steep segment of the low-frequency 
shoulder, the emission from the power-law population (which is more efficient at low frequencies) 
exceeds that of the thermal electrons, while absorption is still mostly dominated by the more 
numerous thermal electrons. As a result, this segment of the spectrum assumes a shape roughly 
described by the rather unusual source function S = Jpi/ath- At still lower frequencies, where the 
emission comes from larger radii in the flow, the contribution of the non-thermal electrons to the 
absorption becomes non-negligible and the spectrum falls more steeply with decreasing v. 

In the region of the thermal peak, we have both jth 3> jpi and ath 3> ctpi, with ath-R ^ 1 
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for frequencies smaller than the peak frequency and ath-R <C 1 for frequencies higher than the 

peak frequency. The spectrum is essentially the same as for a purely thermal model (dashed 
line). (But, note that for sufficiently high r], the synchrotron spectrum assumes an entirely 
non-thermal character and this region too can be dominated by the power-law electrons, with the 
peak luminosity becoming higher and the peak flattened). Beyond the thermal peak, there is no 
self-absorption either by thermal or power-law electrons {a^iR, ap\R ^ 1) and the emission is 
optically thin. Here we find an extended power-law tail dominated entirely by the non-thermal 
population. This segment of the spectrum has the familiar form uL^, ~ jy-(p-3)/2 -^j^jch gives a 
rising spectrum with increasing frequency for p < 3 and a falling spectrum for p > 3. 
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Fig. 3. — Variations in the synchrotron spectrum as a function of (a) the fraction of non-thermal energy content 
T], (b) the power-law index p, (c) the black hole mass in units of 10® Mq, and (d) the mass-accretion rate in units of 
lO-^MEdd- 



Figures 3a-d show how the spectrum of a hybrid model depends on the various parameters. In 
Figure 3a, the energy content of the power-law population r] is varied from 0% (purely thermal) to 
10%; we see that the range and normalization of the low- frequency shoulder increase linearly with 
increasing 77. The normalization of the high-frequency tail also depends linearly on the number of 
electrons in the power-law distribution and thus varies in an obvious way with 77. 

In Figure 3b, we fix rj at 1% and vary the power-law index p. Note that the overall luminosity 
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of both the low-frequency shoulder and the high-frequency tail show considerable dependence on 
p. This is not unexpected. When p is small, there are more particles both at intermediate Lorentz 
factors (log 7 = 1 — 2), which give rise to the low- frequency shoulder, as well as at large 7, which 
emit the high-frequency tail (see §4). Therefore, small values oi p give more emission in both these 
segments of the spectrum. The high-frequency tail has a spectral slope equal to —{p — 3)/2 and 
therefore has a strong dependence on p. 

Two striking qualitative results in Figures 3a and 3b need to be highlighted. First, even 
very small values of the non-thermal energy content ij give rise to significant excess super-thermal 
emission at low frequencies. For instance, even if only 1% of the total electron energy is in 
non-thermal electrons, there can be several orders of magnitude higher luminosity at frequencies 
below the thermal peak. Even more striking is the fact that the spectral shape of the low 
frequency shoulder is nearly independent of rj and p, when both these parameters are independent 
of radius, while the normalization depends on both (see Appendix A for an analytic expression for 
the spectral slope of this segment) . This leads to a degeneracy between rj and p in models of the 
low-frequency emission, so that it is not possible to distinguish between different combinations of 
r] and p, by studying spectral data below the thermal peak only. The degeneracy can be lifted 
with data on the high-frequency tail whose slope has a strong dependence on p. 

Wc finally consider the effect of changing the black hole mass and the mass-accretion rate. 
We find that the relative excess emission due to non-thermal electrons increases rather weakly 
with increasing black hole mass (Figure 3c). This is because the thermal peak moves to lower 
frequencies for higher me (Mahadevan 1997) and at these frequencies the effect of non-thermal 
electrons is more prominent. The non-thermal contribution to both the low and high frequency 
emission increases only weakly with increasing m_3 (Figure 3d). 

3.2. Application to Sgr A* 

Over the last decade, many models have been developed to explain the radio spectrum of 
Sgr A*. This source is believed to be an accreting supermassive black hole at the center of 
our Galaxy. Nearly all the published models invoke synchrotron radiation from relativistic or 
quasi-relativistic electrons. Melia (1992) considered emission by thermal electrons in a spherical 
accretion flow and showed that the resulting cyclo-synchrotron emission is consistent with the 
broad features of the observed spectrum. Narayan, Yi, & Mahadevan (1995) and Narayan 
et al. (1998a) developed ADAF models of Sgr A* which included rotation, viscosity and a 
two-temperature plasma, and obtained similar results, again with a purely thermal distribution of 
electrons. Beckert &: Duschl (1997) and Falcke & Bicrmann (1999) considered non-thermal models, 
while Mahadevan (1998) analyzed a specific hybrid model in which the non-thermal electrons are 
produced by pion decay. This section is a generalization of Mahadevan's work. 

In Figure 4 we apply our hybrid emission model to Sgr A*. We take niQ = 2.5 (Eckart 
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Genzel 1997; Ghez et al. 1998) and adjust m_3 in order to fit approximately the thermal peak. 

The data shown are the same as in Narayan et al. (1998a). Figure 4 shows four models with 
p = 2.5, r] = 0.05%; p = 3.0, r] = 0.2%; p = 3.5, r] = 0.5%; and p = 4.0, r] = 1%. We first note 
that the agreement with data at low frequencies is significantly better with these hybrid models 
(dashed lines) than with a purely thermal model (solid line), as was first shown by Mahadevan 
(1998). 

The two major results pointed out in §3.1 are evident in this figure. First, there is no 
unique solution for the parameters p and rj. The four hybrid models shown in Figure 4 give 
indistinguishable spectra below the thermal peak, rendering it impossible to determine p and 
r] from low-frequency spectral data alone. Second, rj is extremely small in all models. A very 
small fraction of the energy in non-thermal electrons, with rj at most 1%, is sufficient to produce 
all the observed emission at low frequencies. This means that the non-thermal electrons are a 
minor perturbation on the electron population, which is itself a minor perturbation on the more 
dominant ion population. We are therefore consistent when we compute the gas dynamics with a 
purely thermal model and ignore the non-thermal electrons for the dynamics. Further justification 
of this assumption as well as implications of this constraint are discussed in §5. 

We can obtain additional constraints on the parameters rj and p in Sgr A* by studying the 
infrared data. Sgr A* is quiet at infrared wavelengths, with a current upper limit of 10^^ erg 
on the luminosity at 2.2 fim (Eckart & Genzel 1997; Ghez et al. 1998). Since uL^ decreases with 
increasing frequency for p > 3, we cannot use IR data to constrain rj very strongly in these cases. 
On the other hand, if p < 3, electrons in an extended power-law distribution produce significant 
emission in the infrared. For example, Hp = 2.5, we find that the maximum allowed fraction of 
energy in non-thermal electrons is rj = 0.05%. Tighter bounds on the infrared flux will constrain 
the parameters rj and p even more strongly. 

We note that imposing a maximum Lorentz factor 7max on the power-law electron distribution 

also has the effect of suppressing the high-frequency emission. Therefore, for p < 3, we could 
alternatively use the IR data to constrain 7max rather than rj. The 2.2 jim. emission is produced 
predominantly by electrons with log 7 ^ 3 (see §4) placing a maximum Lorentz factor at 
7max ~ 10^ if p < 3 and 77 > 0.05%. 

Finally, we have considered ADAF models with strong outflows following the ideas described 
in Blandford & Begelman (1999), Di Matteo et al. (1999), and Quataert & Narayan (1999). We 
flnd that the constraints on Tj and 7max 

obtained above do not depend strongly on the presence of 

winds. 
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Fig. 4. — Comparison of spectral models for Sgr A* with radio and IR data. The dashed curve shows the spectrum 
when the electrons are purely thermal. The other four curves show spectra from hybrid populations with the following 
combinations of parameters: p = 2.5, r] = 0.05%; p = 3.0, r] = 0.2%; p = 3.5, tj = 0.5%, and p = 4.0, r] = 1%. 



Wc now investigate the effect of a power-law electron population on the size and shape of 
the radio image of an ADAF. As a typical example, we consider an ADAF model of Sgr A* with 
me = 2.5 and m_3 = 0.1 and take a hybrid electron distribution with p = 3.5, rj = 0.5% which 
agrees well with the observed spectrum as shown in Figure 4. Preliminary size measurements 
of Sgr A* are available at two frequencies (Lo et al. 1998, Krichbaum et al. 1998) and shape 
measurements may be possible in the near future. We note that although we use the parameters 
for a specific source, the qualitative results apply equally well to other massive black holes with 
different accretion rates. 

We define the image size as twice the radial distance from the center of the image to the 
point at which the specific intensity falls to half the value of the central intensity, i.e., the FWHM 
of the radio map. Figure 5 shows the predicted FWHM of the image as a function of frequency. 
The dashed curve corresponds to the case when all the electrons are thermal and the solid line 
to the case of a hybrid energy distribution. The two data points for Sgr A* are from Lo et al. 
(1998) at 7 mm and Krichbaum et al. (1998) at 1.4 mm. Note that our radiative transfer code 
is one dimensional and can only handle spherical models, though ADAFs in general are oblate 



3.3. Image sizes and shapes 
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and are likely to appear elliptical in projection (Narayan & Yi 1995a). In the comparison with 
observations, we use only the measured long axes of the images. 




log(Frequency / Hz) 

Fig. 5. — The solid line shows the FWHM image size corresponding to one of the models described in Figure 4, 
with p = 3.5 and q = 0.5%. This is compared with two size measurements of the Galactic Center source Sgr A* 
at 1.4mm (v = 214 GHz) and 7mm {v = 43GHz). We do not show error bars of the 7mm measurement (Lo et al. 
1998) because of the oblateness of the observed image which we have not included in our models. The dashed line 
corresponds to a pure thermal model. 

The effect of power-law electrons is to enlarge the image size at long wavelengths {u < 10^^ Hz), 
where the excess non-thermal emission, which extends outward of the surface of unit optical 
depth, is most prominent. This is a general result which holds true for all model parameters we 
have studied. The other related result is the steepening of the dependence of the image size on the 
frequency at long wavelengths. While for thermal electrons we find 

FWHM/Rsch - ^ (27) 

when we include power law electrons, we find 

FWHM/Rsch ~ ^ (28) 

This is in agreement with the preliminary result obtained by Lo et al. (1998) when they combined 
their measurement of the intrinsic size at 7 mm with the results of Krichbaum et al. (1998) at 
1.4 mm, although the thermal and hybrid slopes are not distinguishable with the current data. 

We also calculate the effective brightness temperature Ti, at each frequency as 
Tb = L^c^/(87r2A;sz^2FWHM^). Figure 6 shows the result. For the thermal model, Tb 
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measures the electron temperature at the photosphere. Since the photosphere moves out 
with decreasing frequency, the brightness temperature falls. For the hybrid model, at lower 
frequencies, the contribution of the non-thermal electrons is large, the spectrum is non-thermal, 
and the resulting brightness temperature increases. This is another clear signature of a 
synchrotron-emitting hybrid electron population. 




log(Frequency / Hz) 

Fig. 6. — The predicted brightness temperature of Sgr A* as a function of frequency. The fiux was taken from 
Figure 4 and the size estimate from Figure 5. The dashed line corresponds to the purely thermal model and the solid 
line to a hybrid model with p = 3.5 and ij = 0.5%. 

We finally study the effect of the non-thermal electrons on the shape of the radio image. 
Figure 7 shows the variation of the brightness temperature Tj, across the source, for the thermal 
and hybrid models of Sgr A* described above. We consider two wavelengths, 3.6 cm and 7 mm. 
The most striking feature of the hybrid case is the limb brightening seen at long wavelengths, a 
feature which is absent in the thermal models. Thus, if the accretion flow contains non-thermal 
electrons, its image would look like a shell rather than a disk. This is again due to the different 
radial dependences of absorption (which is predominantly thermal) and emission (which is mostly 
non-thermal) at the frequencies where the non-thermal shoulder appears in the energy spectrum. 
Because the total absorption falls off more steeply away from the center than the total emissivity, 
the image appears brighter for a range of impact parameters away from the center than at the 
center. In addition, due to the overall increase in the intensity of the emerging radiation in this 
same frequency range, the image looks brighter overall, with the brightness temperature of the 
hybrid case increasing to twice the thermal value in the center of the image at 3.6 cm. The limb 
brightening, along with the enhanced overall brightness at long wavelengths are signatures of a 
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non-thermal population which may be accessible to observations. Unfortunately, in the case of 

Sgr A*, precisely at the wavelengths where the effects arc strongest, interstellar scattering blurs 
the observed image (e.g., Lo et al. 1998). It may be worthwhile to look for these effects in other 
sources where the scattering is less severe. 



Fig. 7. — The predicted image shape of Sgr A* at 7mm and 3.6crn. The dashed line corresponds to the purely 
thermal model and the solid line to the hybrid model with rj = 0.5% and p=3.5. Note that the hybrid model has a 
higher brightness temperature at both wavelengths, and is significantly limb-brightened at 3.6cm. 



4. Electron Energy Distributions and the Shape of Synchrotron Spectra 

The numerical studies presented in §3 show that the presence of a small population of 
power-law electrons causes universal modifications to the spectrum: it introduces a shoulder at 

low frequencies and a power-law tail at high frequencies. We now attempt to associate each of 
these features of the hybrid spectrum with a specific range of Lorentz factors of the non-thermal 
electrons. In doing so, we address three important issues. 

The first is the question of degeneracy: why is the low frequency shoulder in the spectrum 
degenerate to a combination of the power-law index p and the energy content r] of the non-thermal 
electrons? Understanding the source of non-uniqueness is especially important in trying to extract 
the underlying electron distributions from spectral data and determining what we can conclude 
uniquely about these distributions. 
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Second, since introducing even a very small fraction of power-law electrons results in 

significant enhancement of the luminosity at low frequencies, existing radio data can be used to 
constrain the non-thermal energy fraction in these accretion flows. If only a part of the electron 
energy distribution is responsible for the excess emission, data can further constrain these specific 
parts of the particle distribution. 

The third question we address here is related to the absence of significant IR emission from the 
accretion flow around Sgr A* (Eckart &; Genzel 1997), in the frequency range where the extended 
optically-thin non-thermal emission is expected. Absence of emission at these wavelengths provides 
information about the underlying electron distributions by constraining either the power-law index 
p or the maximum Lorcntz factor "jmax electrons can attain in these accretion flows. Determining 
Imax in a- particle distribution may provide a better understanding of particle acceleration and 
cooling processes and timescales in hot accretion flows. 

4.1. Correspondence between electron distributions and photon spectra 

We first investigate the relationship between the Lorentz factor 7 of an electron and the 
dominant frequency of the synchrotron emission it produces in the presence of a dominant thermal 
electron population. At very low frequencies, synchrotron emission comes from large radii in 
the flow where the temperatures are low, the non-thermal emission and absorption are efficient, 
and Opi ~ ath- However, at higher frequencies near and above the low- frequency shoulder in the 
spectrum, ath ^ ctpi, as discussed earlier. For these frequencies, we can thus neglect ap\. We make 
this key simplification only in this section. It simplifies the radiative transfer equation and allows 
us to compute separately the contribution to the intensity from each value of 7 in the non-thermal 
distribution. 

We proceed by writing all quantities as integrals over the power-law electron distribution, 




(29) 




(30) 



and 



(31) 



where nP'(r, 7) is assumed to be a separable function of r and 7 given by 



nf{r, 7) = <(r)<(7) = A^pi(0(P - ^h'"- 



(32) 



(Note that, for this analysis, the non-thermal distribution docs not have to be a power-law, though 
we have assumed this in order to compare this analysis directly to our numerical results). 
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Fig. 8. — (a) The contributions to the hybrid spectrum of individual electron Lorentz factors: 7 = 5, 10, 30, 100, 
and 1000. (b) The range of Lorentz factors contributing to the emission at 10® Hz (low- frequency shoulder), 10^^ 
Hz (thermal peak) and 10^^ Hz (high-frequency tail), (c) A further quantification of (b) showing the minimum and 
maximum Lorentz factors at which the integrand falls to half its maximum value at each frequency. 
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Substituting equations (29)- (32) into the radiative transfer equation and rearranging terms, we 
obtain 

/ [-"^5^ + «th(r)/^(7) - n?\r)j,{j) - jth] n%)d^ = 0. (33) 

For this integral to vanish for any non-thermal electron distribution, the integrand must be 
identically zero, thus giving 

/^^^ = -«th(r)/,(7) + nf{r)M^) + jth- (34) 

We solve this equation for a wide range of Lorentz factors, log 7 = 0.1 — 3.6, by choosing a specific 
radial profile of the non-thermal electron density iVpi(r) corresponding to p = 3.5,1] = 0.5%. We 
study the contribution of each Lorentz factor to the different parts of the spectra. The total 
intensity for the power-law distribution takes the form 

= {p-l)J j-P+'l^{j)d{logj), (35) 

and therefore, in order to study the true contribution of each 7 to the total spectrum, we plot 
the integrand 7~^'''^/7(7) as a function of frequency for each Lorentz factor. The result for 
7 = 5, 10, 30, 100, and 1000 (for p = 3.5) are shown in Figure 8a. The curves give a good 
estimate of the individual contribution of each Lorentz factor to the total spectrum. Most of the 
contribution to the low- frequency shoulder comes from electrons with 7 ~ 30 — 50. Emission from 
electrons with 7 ~ 100 is already lower by an order of magnitude at those frequencies, and the 
emission completely dies off beyond 7 = 100. Contribution to the high-frequency tail, on the other 
hand, starts around 7 ~ 100 and increases with increasing electron energy. 

To get a more quantitative idea of which range of Lorentz factors contributes to the three 
distinct regions of the spectrum, we show in Figure 8b the above integrand as a function of 7 
for 3 frequencies, namely 10^ Hz (in the low- frequency shoulder), 10^^ Hz (in the thermal peak), 
and 10^^ Hz (in the high-frequency tail. We see that the emission at 10^ Hz is primarily from 
electrons with 7 ~ 10^'^; the emission at 10^^ Hz is mostly from 7 ~ 1 (thermal electrons, while 
the emission at 10^^ Hz comes mostly from 7 ^ 

To quantify this effect further, we plot in Figure 8c the minimum and maximum Lorentz 

factors for which the value of the integrand drops to half its maximum value for each frequency. 
This gives us exactly the contributing range of Lorentz factors to the emission at each frequency 
in the low-frequency shoulder and the high-frequency tail. 

We may summarize the results as follows. The low- frequency shoulder is caused by a narrow 
range of electron Lorentz factors, log 7^1 — 2. The narrowness of the range explains why the 
spectra are degenerate to different combinations of p and rj: a power-law distribution with small 
values of p and small rj has nearly the same number of electrons in this narrow range of Lorentz 
factors as one with a larger p and a larger r], thus producing the same emission in the shoulder. 
The high-frequency tail, however, is produced by electrons with a wide range of Lorentz factors 
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(log 7 ^ 2) with higher frequencies coming from higher 7. Intermediate Lorentz factors log 7 ^ 2 
emit predominantly at frequencies around the thermal peak and are overpowered by the thermal 
emission, and therefore do not have an observable feature in the hybrid spectrum. 



5. Energetics of the thermal and non-thermal populations 

In this section, we study the energy flow through the thermal and non-thermal electron 
populations as well as the energy exchange between these two populations through synchrotron 
self-absorption and Coulomb collisions. Starting with the energy equations for the two populations 
and a parametrization of the energy input into each, we calculate the steady-state energy content 
of both as a function of radius assuming that the shapes of the distributions are known a priori. 
Note that if the acceleration mechanism is known, the particle distributions can be calculated 
exactly (see, e.g., Nayakshin & Melia 1998). Here we simply assess the feasibility of a steady-state 
thermal/power-law hybrid distribution. We then extend our discussion of hybrid synchrotron 
spectra to cases where the parameter ry is not constant but is allowed to vary with radius. 



5.1. Energy Equations 

We start by writing the energy equations for the two populations. Neglecting advection 
and diffusion which we estimate to be minor corrections, the energy balance for the power-law 
population reads 

dt 

while the thermal population obeys 



= dntEyisc{r) - / jrix.{ri,p,r)dvdQ. - (36) 



^^'"^ = 5^visc(r) - / ath{B, - h)dvdn + (37) 



at J ie„ 

where -Eth and E^t are the steady state energy contents of the thermal and power-law populations 
respectively, 5 is the fraction of viscous energy that heats the thermal electrons, and (5nt is the 
corresponding fraction injected into the power-law population ((5nt = ^inj*^, where ?7inj is defined in 
§2.2); the time derivatives are set to zero in steady-state. The energy exchange (thermalization) 
timescale due to Coulomb collisions of high energy electrons with the thermal bath is denoted 
by tee- The radiation term in equation ( ^6|) corresponds to the energy loss by the non-thermal 
electrons via optically thin synchrotron emission; because the non-thermal absorption coefficient 
is negligible throughout the flow, the non-thermal electrons do not gain energy by absorption of 



synchrotron radiation. The corresponding term in equation (37), on the other hand, describes the 
heating of the thermal electrons by the local radiation field that is in excess of the blackbody limit. 
When the radiation energy density is a blackbody, there is locally no heating or cooling, but as 
the radiation energy density is above this limit due to emission by the non-thermal electrons, the 
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thermal electrons are heated by absorbing this total emission. Note that this is the same energy 
exchange mechanism invoked in the synchrotron boiler process (Ghisellini, Guilbert, & Svensson 
1988; Ghisellini, Haardt, &: Fabian 1993; Ghisellini, Haardt, & Svensson 1998). The transport of 
energy due to non-local synchrotron self-absorption is likely to be small and have little effect on 
the spectra as the photon mean free path becomes very short very rapidly along a radial path and 
thus the radiation emitted in the optically thick regions of the flow is reabsorbed locally. This is 
due to the very steep radial dependence of the synchrotron absorption coefficient. 

There are several issues we would like to address regarding the relative importance of thermal 
and non-thermal electrons in the flow and the energy exchange between the two populations. First, 
the quantity [5nt/<5](?') determines as a function of radius the relative rate of heating of non-thermal 
electrons compared to thermal electrons. Second, the energy exchange between the populations 
(heating of the thermal electrons by non-thermal electrons) proceeds both via the Coulomb term 
and the absorption of the synchrotron photons emitted by the power-law populations (the second 
term of equation (|37[). Therefore the magnitude of these terms relative to the viscous heating of 
the thermal electrons needs to be assessed. Finally, the importance of the thermalization of the 
power-law distribution due to the Coulomb term also needs to be understood. 



We first estimate the relative magnitude of the terms in equation (pq). The ratio of the 
energy loss of the power-law electrons due to Coulomb collisions (the last term) to the energy 
emitted in synchrotron photons (the second to last term) determines whether Coulomb collisions 
or synchrotron self-absorption is the dominant mechanism by which the non-thermal electrons 
cool. The hybrid spectra presented in this paper have been computed under the assumption that 
the power-law electrons lose energy primarily through synchrotron emission, and this needs to be 
checked for consistency. For a given electron velocity (3 in the power-law distribution, the energy 
exchange rate with the thermal electrons is given by (Nayakshin & Melia 1998) 

'dE\ _ 3 In A 1 

'dtJ^^^ll^P^' ^^^^ 
where InA « 20 is the Coulomb logarithm, is the average thermal Lorentz factor, 
tx = {neCaT)~^ is the Thomson mean-free time and rie the electron density. The rate of energy 
loss of an electron of Lorentz factor 7 emitting synchrotron radiation in a region of magnetic field 
B is (Rybicki & Lightman 1979) 



(39) 

oin-c" 

syn 

where ti/2 is the time for the electron to lose half its energy. Using the analytic expressions for 
the electron density and magnetic field strength given in §2.1 and substituting values for the 
coefficients and parameters appropriate for Sgr A*, we find for an electron of a given 7: 



,^0.1r7-^ (40) 

(£//tl/2jsyn 

where r is the radius in Scwarzschild units as usual. Most of the emission in the low-frequency 
shoulder in the spectrum originates from around r ~ 100 and the radiation is produced by 
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electrons with 7 ~ 50 (cf. Fig. 8c and the Appendix). We thus see that the ratio in equation (40) 
is at most a few percent. This shows that Coulomb collisions play a negligible role in the cooling 
of non-thermal electrons. It also shows that that the heating of thermal electrons by Coulomb 
collisions is unimportant compared to energy exchange via synchrotron emission and absorption. 
We therefore neglect the Coulomb term in the following calculations. 
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Fig. 9. — The ratio of the non-thermal electron heating parameter Snt to the thermal electron heating parameter 
(5 as a function of radius for a model of Sgr A* with p = 3.5, r; = 0.5%. 



The power-law electrons then obey a simple energy balance between the injected energy 
through viscous dissipation and the energy loss through synchrotron emission. The latter simply 
is the total frequency- and angle-integrated synchrotron emissivity of the non-thermal electrons 
because this population is optically thin to synchrotron emission. Therefore, once Sjitir) is 
specified, it is possible to calculate the parameter -q, which describes the steady-state energy 
content of the power-law electrons, as a function of radius. Conversely, if we specify r/(r) as we did 
in the previous sections, it is possible to compute the energy 6nt that would need to be injected 
into the non-thermal electrons as a function of radius. We first study the latter case for t] = 
constant, as this is the assumption we have made in most of this paper. Figure 9 shows 6nt/S as 
a function of radius for one of the hybrid models of Sgr A* with r] = 0.5%, and p = 3.5. The 
figure demonstrates that the energy input into the non-thermal electrons never exceeds 20% of the 
heating of the thermal electrons and in fact does not exceed 10% at those radii that contribute 
to the low-frequency shoulder in the spectrum {lORs < R < lOOOi?^). Thus, at no radius in the 
flow do the power-law electrons become energetically more important than, or even comparable 
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to, the thermal electrons. This may seem like a surprising result considering that the non-thermal 
radiation clearly dominates over the thermal emission in the low-frequency shoulder of the 
spectrum and in the high-frequency tail. However, the synchrotron emissivities of the thermal and 
power-law populations are different by many orders of magnitude at these frequencies. (Optically 
thin synchrotron emission from a power-law population peaks at a much higher intensity and 
at a lower frequency than the thermal emission for the same energy density and magnetic field 
strength.) Therefore, even when the non-thermal population has less energy content than the 
thermal population, and even when a fraction of the non-thermal radiation is absorbed by the 
thermal electrons at frequencies where the flow is optically thick, the escaping radiation can still 
be dominated by the non-thermal emission and the intensity can be much above the blackbody 
limit. This is the situation in the low-frequency shoulder of the spectrum. In the high-frequency 
tail, the only particles that are energetic enough to produce the radiation are the non-thermal 
electrons, and since the flow is optically thin all the radiation escapes freely. Incidentally, the fact 
that the energy input into the non-thermal electrons does not exceed 20% of the heating of the 
thermal electrons (and even this level is reached only over a small range of radius close to the black 
hole) shows that the dynamics of the flow is not affected by the presence of non-thermal electrons. 

We also consider the other mechanism of energy exchange between the two populations, 



namely synchrotron self- absorption by thermal electrons. The second term in equation (36) 
measures the maximum energy that can be radiated by the non-thermal electrons as a function 
of radius. This term also represents an upper limit on the energy that can be transferred to the 
thermal population via self-absorption. This term is bounded from above by ^nt-E'visc according 
to equation (|3^). Since (5nt/5 never exceeds 20% as shown in figure 9, heating of the thermal 
electrons by absorbing synchrotron emission from non-thermal electrons can at most introduce a 
20% correction to the thermal energy equation; the correction is only a few percent at large radius. 
This again justifies neglecting the contribution of the non-thermal electrons to the thermodynamic 
properties of the flow. 

Finally we note that the calculations in the previous sections considered a simple 
parametrization of the non-thermal population via the quantity r], which measures the fraction of 
the electron energy in steady state that is present in non-thermal electrons. Now, rj is a secondary 
quantity whose value depends on the balance between non-thermal heating and cooling. It would 
perhaps be more useful to parametrize the model by specifying the non-thermal heating parameter 
6nt- To this end, we consider models in which the energy injection rate varies as a power-law 
in radius, 5nt r~'', and show results for q between and 1. We pick the normalization of 5nt 
at the inner edge of the flow such that the resulting r/(r) is comparable to the values used in §3 
and is in accord with the results of Figure 9. In Figure 10a, we show r/(r) for p = 3.5 and three 
different choices of q: 0, 1/2, and 3/4. For none of the three cases do we see very large variations 
in the resulting r]{r) versus radius; we flnd that r/(r) rises with radius for q ^ 1/2 and falls for 
larger values of q. Figure 10b shows the spectra corresponding to the three models of (5nt- We 
see that the two universal features of a hybrid spectrum, the low-frequency shoulder and the 
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high-energy tail, are present for all values of q. Thus, super-thermal emission is not an artifact 

of specific assumptions or parameter values of our models but is indeed a robust signature of a 
hybrid population. The slope of the low-frequency shoulder depends very mildly on the particular 
model of 6{r). 




log(R / Rg^j) log(Frequency / Hz) 

Fig. 10. — (a) The steady-state non-thermal energy content 'q{r) corresponding to three models of the non-thermal 
energy injection rate, Snt oc r~^, with g = 0, 1/2, and 1. All the models have p = 3.5. (b) Spectra corresponding to 
the same three models, compared with data on Sgr A*. 



6. Discussion 

In this paper we considered hot accretion flows around supermassive black holes, using the 
ADAF model as a typical example. Such hot flows are expected to occur at low mass-accretion 
rates. We assume that a fraction of the viscous dissipation energy in the accretion flow goes 
into accelerating electrons to a non-thermal power-law distribution. We find that a power-law 
tail of high energy electrons can be sustained in such flows; neither Coulomb collisions nor 
synchrotron self-absorption is able to thermalize the power-law electrons. Assuming that there 
are no other thermalizing mechanisms (e.g., collective plasma modes), we calculate the resulting 
hybrid thermal/non-thermal spectrum from such a plasma. 

The presence of even a population of non-thermal electrons gives rise to two universal and 
prominent features in the synchrotron spectrum: a low-frequency shoulder and a high-frequency 
tail (Fig. 2). These features were identified by Mahadcvan (1998) who considered a specific 
mechanism (via pion decay) for the production of power-law electrons. Even if only a small 
fraction of the total steady-state electron energy is in the non-thermal power-law component, we 
find that there is significant super-thermal emission in the low-frequency shoulder and the high 
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frequency tail of the spectrum. Furthermore, each of these universal features can be associated 
with a specific range of Lorentz factors of the emitting electrons. The low-frequency shoulder is 
emitted by electrons with Lorentz factors log 7 = 1 — 2, while the high-energy tail is emitted by 
electrons with log 7 ^ 2 (Fig. 8). 

Since the power-law electrons cause significant emission at low frequencies, comparing this 
low-frequency shoulder to data can provide stringent constraints on the fraction of the electron 
energy that is present in a non-thermal population. In the case of Sgr A* (Fig. 4), we conclude 
that, in steady-state, electrons even with intermediate Lorentz factors (7 as low as 10) make up 
only a very small fraction of the electron population. Since Coulomb collisions and synchrotron 
self-absorption are ineffective in thermalizing these flows, this means that either the energy 
injected into a non-thermal population, SatE^isc, is small as in the calculations described in §5 or 
there is some other efficient thermalization mechanism which plays a role in these flows. The study 
on Sgr A* needs to be extended to other sources before we can judge how universal the result is. 

Here it is important to emphasize two points regarding the constraints on non-thermal 
electrons obtained from the spectrum. 

First, for the part of the spectrum below the thermal peak, there is a mapping between each 
frequency in the spectrum and the thin shell of radii which predominantly emits at that frequency. 
This is due to the steep radial dependence of the synchrotron absorption that causes a sharp 
transition between optically thick and optically thin regions in the flow. Therefore, emission at 
each frequency comes primarily from a specific radius in the flow and can be used to probe the 
non-thermal content in a shell that extends from that radius out to a somewhat larger radius. 
(Mathematically, radiation emerging at each frequency comes from the entire range from the 
radius where the optical depth is unity out to infinity, but because the total emissivity also drops 
rapidly with radius, the biggest contribution to the intensity integral comes from a relatively 
narrow shell). As a result, the constraints on the non-thermal distribution obtained from each 
frequency of the low-frequency shoulder apply only to electrons within this shell. Inside the 
radius corresponding to optical depth unity, the emission is heavily self-absorbed and the electron 
distribution there is inaccessible to observations. 

Second, the constraints on the non-thermal energy content which we have obtained for Sgr A* 
would not be invalidated if the observed low-frequency shoulder originated in a separate, extended 
region such as a jet or an outflow exterior to the accretion flow. In that case, the observed fluxes 
may be treated as upper limits to the emission from the inner accretion flow. This would imply 
that the non-thermal emission from the accretion flow is even smaller than in our models and the 
resulting constraints on the fraction of non-thermal electrons would become only tighter. 

When the non-thermal electron distribution is a power law, three parameters of the model 
and three corresponding pieces of physics can be extracted by comparing hybrid emission models 
to spectral data. The first parameter is rj, or equivalently ^nt (cf.§5), which measures the fraction 
of the dissipated energy that goes into the non-thermal electrons. The second parameter is p, the 
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slope of the power-law energy distribution (eq. 4). This can be uniquely determined by measuring 
the spectral slope of the high-frequency tail in the spectrum (provided there is no competing 
emission via Compton scattering at these wavelengths, cf. Fig. 1). The value of p may provide an 
indication of the mechanism by which the electrons in a hot accretion flow are accelerated into 
a power-law tail. The third parameter, which is more significant for distributions with p < 3, is 
7max, the maximum Lorentz factor to which electrons are accelerated. If this parameter can be 
determined by an IR cutoff in the spectrum, it would provide information on the ratio of the 
acceleration timescale to the synchrotron cooling timescale. 

Due to the degeneracy of the low- frequency shoulder to the parameters rj and p (Fig. 4) , this 
segment of the spectrum by itself is not sufficient to constrain the two parameters. However, j] 
and p may be determined uniquely if we have information on both the low-frequency shoulder and 
the high-frequency tail. For this to work, there should be no contamination from a jet or outflow 
to the observed radiation in the low-frequency shoulder or from Compton scattering in the high 
frequency tail. 

We showed in this paper that the presence of a non-thermal population of electrons also has 
measurable effects on the shape of the image of the source, and the size of the accretion flow, as a 
function of frequency (Figs. 5-7). Because the frequency dependence of the image size is strongly 
affected by the fraction of non-thermal electrons in the flow, measuring the size of accretion flows 
at multiple wavelengths can provide a quantitative constraint on the non-thermal electron energy 
content. Moreover, the brightness temperature of the source as a whole, as well as the variation 
of the brightness temperature across the source, at long wavelengths behave differently for a 
hybrid plasma compared to a purely thermal population of electrons. Finally, at the frequencies 
corresponding to the low-frequency shoulder, we show that there would be limb brightening of 
the image. All these image-related signatures may help to identify the presence of a non-thermal 
population of electrons, and the parameters of the corresponding energy distribution. 

We thank Eliot Quataert, Anthony Aguirre, Tiziana Di Matteo, and James Cho for many 
useful discussions and comments on the manuscript. We also thank S. Nayakshin, the referee, for 
many useful comments and suggestions that improved the presentation of the paper. D. P. was 
supported by a post-doctoral fellowship of the Smithsonian Institution. This work was supported 
in part by NSF Grant AST 9820686. 

A. Analytic Approximations 

In this appendix, we provide analytic formulae for the slope of the spectrum below and above 

the thermal peak and the relative importance of non-thermal electrons as a function of mass, 
accretion rate, power law index p, and non-thermal energy content rj. For convenience, we use the 
self-similar solution of ADAFs to obtain these analytic estimates. 
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A.l. Self-Similar Advection-Dominated Flows 

We begin by presenting the self-similar ADAF solution developed by Narayan & Yi (1994, 
1995b). The self-similar solution describes the local properties of the accretion flow as a function 
of the black hole mass M, the mass accretion rate M, the radius R, the viscosity parameter a, the 
ratio of gas pressure to the total pressure (3, and the fraction of viscously dissipated energy that is 
advected inwards /. 

In terms of the scalings introduced in §2.1, the height-averaged electron number density ng, 
the magnetic field strength B, and the dimensionless proton temperature 9p = kTp/mp(? of the 
accretion flow are: 

He = TiimQ^msr^^^'^ g cm~^, (Al) 
B = hx-ml'I'^mllr-^l^ G, (A2) 

and 

6*^ = 0.18/3r-\ (A3) 

where 

m = 2.0 X 10^°a-^cr^C3 (A4) 

and 

6i = 2.07 X 10^a-i/2(i _ 0^c[^l\}l\ (A5) 

The coefficients c\ and C3 are defined in Narayan & Yi (1995b) and arc related to the adiabatic 
index of the gas and the fraction of advected energy /. For the cases of interest here, c\ ~ 0.5 and 
C3 ~ 0.3. The remaining parameters were specified in §2.1. 



A. 2. Fundamental synchrotron quantities in self-similar ADAFs 

Wc now use the self-similar solution to obtain expressions for the cyclotron frequency and the 
harmonic number for synchrotron emission. Since we focus our applications to galactic nuclei and 
very large mass black holes where the synchrotron frequencies of interest are in the radio regime, 
we scale the frequency as i^io = z^/10^°. The magnetic field strength is then 

B = 4.85 X lO^mg ^/^mi/g r"^/^ (A6) 

and thus fundamental cyclotron frequency defined in §2.3 becomes: 

z/fc = 1.4 X lOiim" ^/^mVf r-^/^ (A7) 

The ion temperature retains its virial value throughout the flow and is given, in dimensionless 
units Qi = by: 

Oi = 0.61/3c3r-^ = 0.2r-^ (A8) 
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when we set the parameter values to those used in the simulations. The behaviour of the electron 
temperature on the other hand is more complicated and depends on the detailed balance of heating 
and cooling at each radius. Qualitatively, electron temperature starts out virial and equal to the 
ion temperature at large radii (r = 10^), but increases less steeply than the ion temperature at 
smaller radii due to cooling by mainly synchrotron emission. Because it is not possible to derive 
the exact electron temperature as a function of radius analytically, we provide instead numerical 
fits to the average radial dependence of electron temperature 

Oe = 1.65 X r-° ^ (A9) 

Since synchrotron absorption coefficient is a steep function of radius, we can assume that 

most of the optically thick emission comes from a narrow range of radii around the radius with 
optical depth equal to unity. Due to the complicated nature of the expression for the thermal 
synchrotron absorption which does not allow analytic integration, we numerically determine the 
dependence of this r = 1 radius, rti on mass, mass accretion rate and frequency. The best fit gives 

m = 2.5 X lO^mg ^^''mi5z/{o°'^ (AlO) 
Substituting Vti into 1^5 and 9e we get 

V, = 1.4 X lQ^m,"''m'!l\%' (AH) 

and 

ee = 0.06 X mg-^^mI^-V?^6 (A12) 
Finally, we compute the harmonic number defined in §2.3: 

XM = ^ (A13) 

XM = 1-32 X lO^e^'^^-a^^ro"'^^ (A14) 

which provides all the necessary synchrotron expressions as a function of black hole mass, mass 
accretion rate and frequency. 



A. 3. Where does the non-thermal emission dominate: Normalization 

We can now derive analytic estimates for the contribution of the non-thermal emission to the 
low-frequency shoulder and the high-frequency tail. In the optically thick shoulder, self-absorption 
of the synchrotron emission is done predominantly by the thermal electrons, and we find that 
absorption due to non-thermal electrons are negligible down to a frequency of ~ 10^ as discussed 
above. Thus the source function will take the form 



(A15) 



-30- 



where _7th is synchrotron emissivity per unit volume given in §2.3. 

For the analytic approximations, we will use a simplified form of the thermal emissivity. We 
first take the fully relativistic limit derived first by Pacholczyk, which is equivalent to setting a, b 
and c equal to 1 in M{xm)- We further neglect the second and third terms of the sum altogether 
by setting b and c equal to zero as these terms provide only a very small correction unimportant 
for the present purposes. Thus, in simplified form, the ratio of the two emissivities is given by 

^ = r?C;,a(^)i^2(l/^)^(-)-(^+')/':ri/^xp(1.9x]^/=^) (A16) 

Jth n 

where 

= ^ (A17) 

and all other quantities are as defined as before. 

The non-thermal emission dominates over the thermal emission in the spectrum when 
ipi/ith > 1- This ratio can be easily computed using the expressions for XM,de, smd evaluated 
at the T = 1 radius given above. 

The power law emission beyond the thermal peak is optically thin and hence the excess in this 
region depends sensitively on the power-law index as expected and is simply given by the optically 
thin synchrotron emission from the power-law electrons. The amount of emission is determined 
simply by jpi. 

We finally note that, for the small values of ij considered here, thermal emission dominates 
at the thermal peak. The strong dependence of the position and the normalization of the thermal 
peak on m and rh is discussed by Mahadevan (1997). 



A.4. Spectral Slope 

The shape of the low-frequency shoulder is dominated by the emission of the power-law 
population and absorption by the thermal electrons. In order to determine the spectral slope of 
this segment, we therefore first calculate the hybrid source function 

Su = Jpi/ath (A18) 

at the surface of unit optical depth. To convert this into a luminosity, we multiply the source 
function by the area of the r = 1 surface. If we write the shoulder of the spectrum as vLj^ oc v^, 
the spectral slope s is then given by 

s = 2.5 - p/8 - 7.2m%iy^Q^-^^ + Gm^^-^^m^^u^Q^-^^ (A19) 

which gives the canonical value of ~ 4/3 for the mass and mass-accretion rate of Sgr A* and 
scales (weakly) according to the expression above for other masses and accretion rates. One can 
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immediately see that the spectral slope depends very weakly on p, only as p/8, thus demonstrating 
again the independence of the spectral slope of the low-frequency shoulder on the shape of the 
electron energy distribution. 

The slope of the optically thin power law emission is a well-known result, 

s = -(p-3)/2 (A20) 

which in vL^^ for p < 3, flat for p=3 and falls oflF for p > 3. It is independent of the black hole 
mass and the accretion rate. 
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